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-- ■ ABSTRACT 
O 
O 

' We re-examine the characteristic polarization angle sweep of rotation-powered pul- 

sars and calculate the expected deviations from this sweep caused by aberrational effects 

■ and by polar-cap current flow. We find that in addition to the previously known phase 
shift of the entire sweep by A<I> = —Ar/Ri, aberration shifts the polarization angle 

, itself by A^' = — (10/3) (r/ii^) cos a. Similarly, current flow above the polar cap shifts 

the polarization sweep by A^' = (10/3)(r/i?L)(J/JGrj) cosa, potentially providing a 
\ method of directly measuring the magnitude of the current. The competition between 

these two effects produces a potentially observable signature in the polarization angle 
sweep. Although these effects may appear similar to orthogonal mode shifts, they are 

OO ' 

■ an independent phenomenon with distinct observational characteristics. 

o ■ 
o 

::a 

p i' 1. Introduction 

o 

\-{ I The characteristic S-curve polarization sweeps exhibited by pulsar emission have spawned 

^ ' a large body of research, both theoretical and observational. Ever since the highly-successful 

rotating vector model of Radhakrishnan & Cooke (1969), hereafter RC, fitting pulsar polarization 
profiles has been a standard method of constraining the underlying pulsar geometry. Some pulsars, 
however, have polarization behavior which does not fit the standard S-curve. The millisecond 
, pulsars, for example, seem to have noisy and on average flatter polarization sweeps than normal 

(Xilouris, et al. 1998). In attempts to understand these deviations, several relativistic and plasma 
effects have been proposed which would modify the sweep through plasma propagation effects 
(Barnard & Arons 1986; McKinnon 1997), through aberration of the beaming direction from strict 
parallelism (Blaskiewicz, Cordes, &: Wasserman 1991), or through multiple interacting orthogonal 
modes (McKinnon &: Stinebring 1998). 
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Such effects become progressively more important as tfie radius of emission (or the polarization- 
limiting radius) approaches the light cylinder distance, Rl = c/Q. r/R^ is very small in the polar 
cap emission model for pulsars (r ~ R*), a model consistent with most observational constraints 
on r (Rankin 1990). Since Rl is significantly smaller for millisecond pulsars, due to their faster 
rotation, we expect these effects to be proportionally larger. 

However, most of the departures from the pure RC model previously considered have depended 
on effects of second-order or higher in r/Ri, such as the field sweepback perturbation of Barnard 
(1986). The only first-order effect which has been considered is the aberrational shift due to the 
co-rotational velocity, examined by Blaskiewicz, Cordes, Sz Wasserman (1991), hereafter BCW. In 
this paper, we re-examine that effect and analyze an additional first-order effect due to current flow 
above the polar cap. 

This is a phcnomcnological study, guided by the success of the rotating vector model. That 
model contains three chief assumptions: that the observed radiation is beamed along the direction 
of the magnetic field lines, that the polarization of the radiation is at a fixed angle to the radius of 
curvature of the field lines, as is the case for vacuum curvature radiation, and that the underlying 
magnetic field is a dipole. To go beyond this, we relax some of these assumptions in physically- 
motivated ways. Keeping the assumption that the polarization is related to the curvature of the 
field lines, we include aberrational change in the beaming direction and the effects of currents on 
the underlying dipole field. 

After taking these effects into account, we find that aberration delays the phase of the polar- 
ization sweep with respect to the center of emission and, in a plot of polarization angle vs. rotation 
phase, shifts the entire sweep downwards, while current flow shifts the entire sweep upwards. If 
both effects are present, they may lead to sharp jumps in the polarization angle near the edges of 
current-carrying zones. 

2. Polarization Perturbations 

Since we assume that the polarization is related to the acceleration of the emitting particles 
and that the emission is beamed along the direction of motion, the particle velocity field is the core 
quantity examined in this paper. Particles follow the magnetic field, so any perturbations to that 
field are refiected in the velocity. 

We apply the same general method to several perturbations. Starting with a simple dipole 
field, we add a perturbation to find the new velocity field of the emitting particles. From that field, 
we can calculate the radius-of-curvature vector of the particles at each point in space. Then, for 
any pulsar phase, we find the position where the velocity field at a given radius r points towards 
the observer, evaluate the radius of curvature vector at that point, and project it onto the sky to 
find the observed polarization angle. 
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To obtain analytic expressions, these calculations are done to first order in t/Rl- This naturally 
gives the difference /\ip between the normal polarization angle sweep, 'ip{^) and the perturbed sweep, 
= ^($) + AV'($). More details of this procedure are given in Appendix A. 

Two angles define the basic geometry of a pulsar: the angle a between the rotation axis (fi) 
and the magnetic moment (//), and the inclination angle i between the rotation axis and the line of 
sight to the observer {(,), see Figure 1. The angle of closest approach, the minimum angle between 
/i and f2, is then j3 = i — a. 

These calculations are done in magnetic-centered coordinates, defined so that the z axis (6 = 0) 
is along the magnetic moment fx, the rotation axis ft is fixed in the x — z {(j) = 0) plane, and the 
vector to the observer £ is also in the x — z plane at t = 0. 

We have also numerically calculated the expected polarization angle sweeps for each pertur- 
bation examined. Given a modified velocity field, we directly solved for the new emission point, 
computed the radius of curvature vector, and compared it to the rotation axis to find a polarization 
angle. These numerical results confirmed the analytic perturbation theory. 



2.1. Offset dipoles 

One suggested perturbation to the magnetic field is an offset dipole, where the center of the 
dipolar field is offset from the center of the star, e.g. Chen & Ruderman (1993). This has negligible 
effect on the polarization sweep. In an offset dipole, the magnetic origin revolves about the center 
of the neutron star with a radius of 5 < i?*. In turn, in the frame of the magnetic field, the 
observer appears to move in a circle of identical radius. The electrodynamics in these field-centered 
coordinates is the same as in the star-centered case, so the only change in the observed emission is 
that due to the effective motion of the observer. This motion and the resultant changes are tiny, 
on the order of <C 1, where d is the distance to the observer. 



2.2. Polar Field Aligned Current 

Most models of the pulsar polar cap include a current of charged particles streaming along the 
field lines, a current approximately equal to the Goldreich-Julian density moving at the speed of 
light. This current induces a magnetic field of Bi = {r / Rl){J / Jgj)B cos a e^; this field is purely 
toroidal with respect to the magnetic axis. 

From Appendix D, the resulting polarization angle shift is 

Atp = — - — cos all — — sm Uq radians (Ij 

3 Rl Jgj V 40 J 

where is the magnetic colatitude of the (0*'*-order) emission point, see Appendix B. At the 




Fig. 1. — Basic Angles, a is the angle between the rotation axis, 11, and the magnetic moment, /it. 
i is the angle between the observer LOS, ^, and the rotation axis. 
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center of the pulse, has its minimum value, ~ (2/3)|a — i\ = (2/3)|/3|. At the boundary of the 
last closed field line, reaches its maximum, ~ (2/3) y^r/i^^. For r/Ri small, the range of 9q 
is also small, and the perturbation is mostly flat across the observed pulse. 

In terms of the orbital phase, the perturbation is 

. , 10 r J /27 5 ^\ 25 r J . „ . . 2^ 

Aw « — — COS a 1 — -cosp -p:— — sm2Q;smzsm lit. (2) 

^ 3RlJgj V32 32 192 Rl Jgj 

To get this form, wc used the approximation formulae of Appendix C, assumed that was small 
enough that cosilt ~ 1 — (1/2) sin^ ilt, and simplified the constants by taking 28/45 ~ 5/8, which 
is good to 0.4%. Here, t = corresponds to the center of the pulse profile, where the slope of the 
polarization angle curve is greatest. 

In terms of observablcs, r/R^ ~ 1.2° PQ~-|^rioo, where rioo is the emission radius in units of 
100 km and Pq.i is the pulse period in units of 0.1 seconds, so the magnitude of this shift is 
4.0°Po:i(<^/'>^Gj)rioo. 



2.3. Aberration 

As discussed by Blaskiewicz, Cordes, Sz Wasserman (1991), the rotation of the magnetosphere 
itself affects the velocities of particles moving along field lines. Inside the light cylinder, the field 
lines are stationary in the frame co-rotating with the star. However, as particles move along these 
lines, their paths curve both because of the curvature of the field line itself and because the frame 
of those field lines rotates with the star. Even if the field lines were straight, the paths of the 
particles moving along them would be helices in the inertial frame. This extra curvature generates 
a perturbation of the instantaneous particle directions equal to the co-rotation velocity, fl x r. 

To lowest order, this perturbation shifts the polarization angle by 

A^ = -2jl """7^° +0(jl). (3) 
^ Rl sm0o ^Rl' ^ ' 



Since sin is generally tiny, on the order of the polar cap size, = \/r/RL, this perturbation is 
of order ^JrjRL^ a half-order lower than the constant shift produced by the current. 

This shift has the same form as a phase shift of the entire polarization sweep by 

A* = -3^ (4) 

with respect to the geometric center of the pulse phase, as was found previously by BCW. Since 
aberration advances the emission itself by A$ = r/RL due to simple beaming, the total phase shift 
with respect to the emission is A$ = —At/Rl- 

Prom numerical investigations, representing this effect as a phase shift closely matches the 
perturbation. The analytic method presented here only gives the first order correction (in t/Rl) 
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to the sweep; higher-order terms quickly become important, since aberration adds appreciable 
curvature to central field lines which formerly had none. Representing the change as a phase shift 
better reflects these higher-order terms, while giving the same first-order result. 

If Figure 3, we show the numerically-calculated polarization sweep for r/R^ = 0.1; this gives 
a good idea of the character of the perturbation. As shown in Figure 4, when r/Ri = 0.01, the 
analytically predicted first-order shift is nearly identical to the numerical, while at v/Rl = 0.1 there 
is a noticeable difference between the two, which is mostly corrected by casting the perturbation 
as a phase shift. 

Once this phase shift is removed, we are left with a remaining perturbation of 



A* = — -— cosa 1-— sin 6*0 + — — sina sin6'o cos 00 + 0(sin'' 6'o). (5) 
6 Kl \ iO / 12 Rl 

See Appendix G for details. In terms of the pulsar phase, this is approximately 

10 r /3 5 \ 47 r 23 r 

A* Ri — -— -— cos a I - + - cos /? I -I sin asm. 13 — -— sin 2a! sin i siv? VLt (6) 

3 i?L V8 8 J 18 Rl IURl ^ ' 

using the same approximations as for Eq. (2) . 

The constant portion of this polarization angle shift cancels the constant portion of the per- 
turbation due to polar current flow when J/Jqj = 1- However, if J/Jqj 7^ 1, there is a nonzero 
shift, providing a diagnostic of the magnitude of current flow. The remaining linear and quadratic 
terms produce slight changes to the shape of the polarization angle sweep. 



2.4. Relativistic Sweep-back 

The actual magnetic field of a pulsar will differ from the ideal dipole due to the relativistic 
sweep-back of field lines, caused by the induced displacement current. The full vacuum fields are 
given in Deutsch (1955) and most recently restated by Melatos (1997). Admittedly, the region 
surrounding a pulsar is almost certainly not a vacuum, but the morphology of the vacuum field 
should roughly reflect the structure of the magnetic field with conduction currents. 

The difference between the Deutsch fields and the vacuum dipole is second order in t/Rl, 
however, so the first order effects considered here dominate at low altitudes {t/Rl < 0.2). 

In the non-vacuum case, the corresponding perturbation is the magnetic field due to the 
rotation of the Goldreich-Julian density. This field component satisfies V x B^ot = 47r(0 x r/c)r]Gj- 
Since both fi x r/c and rr]Gj/BQ are manifestly first-order in v/Rl, Bj-ot/Bo is second order and 
therefore has been neglected. 
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Fig. 2. — Perturbation to polarization angle, A^, due to a Goldreich-Julian current as a function 
of longitude. All angles are in degrees. The dotted line is the unperturbed sweep, while the solid 
line is the perturbed sweep. t/Rl = 0.1, a = 30°, and i = 35°. 
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Fig. 3. — Perturbation to polarization angle, A^, due to co-rotation effects as a function of longi- 
tude. All angles are in degrees. The dotted line is the unperturbed sweep, while the solid line is 
the perturbed sweep. r/Ri = 0.1, a = 30°, and i = 35°. 
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Fig. 4. — Analytic polarization angle shift due to aberration compared to the numerical ioY r / Rl = 
0.01 and t/Rl = 0.1. a = 30°, and i = 35°. The solid line is the analytic perturbation theory, 
the dot-dash line (for t/Rl = 0.1) is the phase shift derived from the perturbation theory, and the 
dotted line is the numerical result. Casting the perturbation as a phase shift better matches the 
numerical result for large t/Rl. 
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3. Return Currents 

While the perturbations described in the previous section are general effects, this section applies 
those effects to a specific model of the polar cap. To maintain charge balance, the current flowing 
out from the polar cap must be balanced by a return current elsewhere in the system. Here, we 
assume this current flows along the boundary layer of last-closed field lines, surrounding the polar 
cone (Goldreich and Julian 1969). 

Outside this auroral sheath, the net current flow through a loop surrounding the pole is zero, 
so there is no magnetic field perturbation due to the current. Only the aberrational shift remains. 
The field on the open field lines stays the same, containing perturbations from both current and 
aberration, while within the layer itself the current perturbation is quickly eliminated. For more 
detailed analytic models, see Appendix E. 

If this return current layer is illuminated by emission, it will produce a characteristic signature 
in the polarization sweep. Figure 5. The sharp transitions are caused by the current perturbation 
turning on and off as the line of sight passes through the return current layer. 

Although the return current layer is not normally thought of as a site of emission, it may be 
visible either through direct emission, by some form of two-stream instability, or through refraction 
or scattering of radiation emitted within the cone. Since we are more interested in the consequences 
of current flow above polar caps, we leave the precise mechanism of the emission open. If the return 
current is composed of outflowing positrons or ions, the beaming properties would be much the 
same as those of the primary beam, while if the return current interpenetrates the outflow, even 
the identical mechansism would reflect the presence of the return current. Transitions like the one 
in Figure 5 would be strong evidence for emission in this region or radiation transfer through it. 

In addition to these specific effects, a more general conclusion is that inhomogeneities on the 
order of Jqj in the current flow above pulsar polar caps should produce perturbations in the 
polarization angle sweep on the order oi v/Rl- 

4. Discussion 

Given a single polarization angle sweep, at a single frequency, observers can hope to see the 

effects of any current inhomogeneities, small deviations from the standard shape caused by higher- 
order terms, and the phase shift between the center of the polarization angle sweep and the pulse 
itself. The constant offsets caused by current flow would be detectable only through comparison of 
multiple frequencies, which would also help confirm any of the other effects. 

From a more qualitative point of view, these concepts may be used to help understand certain 
difficult objects. For example, the potentially large jumps in the polarization angle when the line 
of sight crosses a return current layer could explain some of the shifts which are normally ascribed 
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Fig. 5. — Polarization sweep for a = 30°, i = 35°, and v/Rl = 0.1, with a return current layer at 
$ = ±25° with a thickness set by A = 50, showing the resultant jumps in the polarization angle 
sweep. 
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to orthogonal modes. Clearly, this result does not explain orthogonal modes in general, as can be 
seen from the single-pulse studies of XXX, but these shifts could masquerade as orthogonal mode 
transitions in individual pulsars. 

For example, in J0437-0415, the magnitude of the "orthogonal" transitions changes with fre- 
quency, ranging from 90° at high frequency to 70° and 40° at lower frequencies, suggesting that 
either these are orthogonal modes combined with propagation effects or, perhaps, that the 90° shift 
is a coincidence and these shifts are instead due to current features in the magnctosphere. 

In addition to polarization jumps, J0437-0415 displays a polarization phase lag in which higher 
frequencies lag lower (Navarro 1997). Combined with the increasing transition amplitude with 
frequency, this suggests that the higher frequencies arise from higher in the magnetosphere than the 
lower frequencies, contrary to the normal radius-to-frequency map. Since most millisecond pulsars 
do not display the pulse widening with decreased frequency that is one of the best arguments for 
radius-to-frequency mapping in normal pulsars, it is certainly conceivable that the structure of the 
emission region is different. In any case, the small range of phase shifts seen argues for a small 
emission region, on the order of a few kilometers. 

5. Summary 

In this paper, we have shown how two physical effects, rotational aberration and current flow, 
affect the polarization angle of observed radiation. The strength of these effects is proportional to 
the height of the emission region. 

In brief, current flow tends to shift the entire polarization angle sweep upwards in a plot 
of polarization angle vs. phase, while aberration tends both to shift the sweep downwards and to 
phase shift it. The largest of these is the aberrational phase shift between the center of the emission 
profile and the center of the polarization sweep. 

The Goldreich-Julian current and aberration produce equal and opposite constant offsets to 
the polarization angle sweep, effectively cancelling each other. Any deviations in the current from 
the fiducial Goldreich-Julian value will break this cancellation, giving a potential indicator of the 
magnitude of current fiow above the polar cap. 

Inhomogeneities in the current fiow structure above the polar cap should be revealed as per- 
turbations, of order t/Rl, in the polarization sweep. The specific case of an auroral return current 
layer would appear as sharp transitions along the edge of the observed pulse, potentially correlated 
with conal emission. These jumps may appear similar to orthogonal modes, but should have an 
amplitude which changes with frequency, if frequency is related to emission radiiis. In the presence 
of current sheets, single-pulse polarimetry should show a shift in both orthogonal modes, provided 
that multiple altitudes do not contribute to the observed radiation at one frequency. This is quite 
different from the existing single-pulse studies (Gil &: Lyne 1995; Gangadhara 1997), which show a 
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simple shift in the relative amplitude of two existing modes. 

We would like to thank Don Backer and Andrea Sonier for illuminating discussions about the 
observability of these effects and Anatoly Spitkovsky for helpful animations of swept-back field 
lines. 



From classical differential geometry, the tangent vector is the unit vector along B, t = B/|i?|, 
while the normal vector n is the unit vector along the radius of curvature. For highly relativistic 
flow, radiation is beamed along the direction of the particle velocity, so at any instant the observed 
radiation comes from the spot in the magnetosphere where the tangent vector i points in the 
direction of the observer. 

For normal vacuum radiation, the polarization vector must lie in the plane perpendicular to 
t; for the particular case of curvature radiation, the polarization is almost entirely along n. For 

plasma processes, this restriction is relaxed, but in the absence of other directional perturbations, 
such as density gradients not along n, we expect that the polarization angle would be related to 
the direction of curvature of the field lines. Density gradients along n, which we expect if polar cap 
pair creation is the main source of the plasma, would likewise link the polarization to the magnetic 
geometry. 

To trace the orientation of the field, we use the binormal vector, defined by b = t x n. This 
is equivalent to tracing the normal, but is computationally more convenient, since the binormal of 
a dipole field is simply a constant vector along e^. 

The observed polarization angle is then set by finding the angle between the normal vector 
and the projected rotation axis of the pulsar. Since the tangent points towards the observer, both 
the normal and the binormal are automatically in the plane of the sky, but the rotation axis must 
be projected, giving fip = A — (A • 

If the angle between the rotation axis and the binormal is then cos t/j' = b ■ ftp and 
simp' = t ■ (b X Ctp). Since the angle made by the normal vector is 90° different from that made by 
the binormal, the true polarization angle is given by 



A. 



Geometric Method 



b-{txU) 

where we have expressed the result in terms of Cl, rather than ft. 



tan lb = — cot 



(Al) 



For a standard dipole field, this gives the normal polarization curve. 



tan'0 = 



sin a sin 



(A2) 



cos a sin z — sin a cos i cos ' 
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c.f. Radhakrishnan &; Cooke (1969). 

In general, we want to know the change in the polarization angle due to a small perturbation 
to the background magnetic field, B = Bq + Bi. This change is due to shifts in both the emission 
point and the binormal itself. 

To first order, the new tangent vector t will be the sum of the initial tangent vector, to, 
and a perturbation ti perpendicular to to- (Any component along to is simply absorbed into the 
normalization.) 

The new curvature vector is 

K = (t • V)t = (to • V)to + (to • V)ti + (ti • V)to = Ko + Ki (A3) 

where k = n/p, kq = no/po = (to • V)to, and 

Ki = (to • V)ti + (ti • V)to (A4) 

When normalized, the component of ki along no is absorbed into the normalization, giving a 
perturbed normal vector of n = no + ni , where 

ni = po{ki - (no • Ki)ho). (A5) 



The perturbed binormal is then b = bo + bi, where 

bi = ti X no + to X ni. (A6) 



The second contribution to the polarization angle change comes from the change in the emission 
point. If the unperturbed field points towards the observer at ro, the perturbed does so at ro + ri, 
where ri can be found by solving 

to(ro) + Vto-ri + ti(ro) =1 (A7) 



or 

ri = -[Vto]-^ •ti(ro) (A8) 

This expression is under determined, as there is typically an entire line of points where the field 
points towards the observer. 

In this paper, we assume that emission at a given frequency arises at (or decouples from the 
magnetic field at) a fixed radius, so we constrain the radius to remain constant. This "lateral" 
motion shifts the emission site to a different field line from the unperturbed one, and, if transverse 
gradients are steep, with a size scale shorter than roughly r'^sinO/RL, this will likely affect the 
emission properties. Near the boundaries of the polar cap, such transverse gradients may become 
important, requiring a more careful treatment. 
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This changes the observed binormal by an additional Sh = V6o " ri , giving a total change of 
Ab = bi + Sh. Since this change in the binormal is perpendicular to the unperturbed binormal, 
the magnitude of the change in polarization angle is simply \Aip\ = |Ab|. 

Getting the sign correct requires comparing the change in the binormal to the direction of the 
projected rotation axis, or 

to • (Ab X Qp) 



bo ■ Qp 



(A9) 



Since all of these steps are linear in the perturbation, multiple perturbations may be consid- 
ered separately, then simply added. In later sections, we consider in turn two different possible 
perturbations to the basic dipole field. 



B. Unperturbed Field 

In (r, 9, (p) polar coordinates centered on the magnetic axis, the unperturbed dipole field is 

2/Lt cos 6 iJ, sin 9 



Bn = 







which gives a tangent vector field of 



. _ Bo _ 1 



cos 9, - sm9, 



where 



AT, = 1 - - sin^ 



1/2 



The curvature vector is then 



-• 3 sin 9 ,^ 1 9 ^, 
Kq = (to • V)to = 73774 (1 - o 9) 



1 



sin 6*, cos 6*, 



giving a normal vector of 



no 



Ko 1 



— - sinO, cos 9, 
2 



\ko\ Ni 

The binormal is purely in the ^direction, 

5o = to X no = [0, 0, 1]. 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



In order to relate this binormal field to the actual polarization angle, we need to know the 
emission point, which we denote ro = (ro, Oq, (po), where the magnetic field points towards the 
observer at {9obs, ^obs)- 
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Since the dipole field contains no 0-component, we must have = 0o6s- Solving for Oq then 
gives 



If 9obs is small, as is usually the case for an observer looking down the cone of open field lines, this 
simplifies to Oq ^ {2/3)9obs- 

If wc take these results and evaluate expression (Al), we obtain the standard polarization 
angle sweep, (A2). 



Evaluating the polarization effects of various perturbations naturally leads to expressions which 
depend on the coordinates of the emission point, usually in the frame co-rotating with the star. 
Since no actual observer is sitting in that co-rotating frame, we have to rewrite the results in 
terms of observables, such as the rotational phase (Jit) and the physical inclination angles {a, i). 
Admittedly, the inclination angles are not directly observable, but they set the geometry of the 
pulsar. 

Unfortunately, this conversion is a significant source of error in the analytic perturbation 
theory. Straightforward power-series expansions may be used to cast the raw analytical expressions 
into more useful forms, but doing so limits the range of validity to a small portion of the pulse 
phase near the pulse center. 

Here, wc first give the exact formulae relating the pulse phase to the position of the observer 
in these coordinates. These, coupled with the solution for the 0-order emission colatitude (B7), 
give a precise definition for Oq and 0o at any pulse phase. We then give several more tractable, but 
more limited, approximations. 

The vector to the observer is 




(B7) 



C. Approximations 



e = 



cos a sin i cos Qt — sin a cos i 

— sin i sin 0,t 
sin a sin i cos fit -|- cos a cos i 



(CI) 



in Cartesian magnetic coordinates. In pola] 



ir, the observer is at {9obsi4>ohs) given by 




cos Oobs = COS a cos i -I- sin a sin i cos Clt 



(C2) 
(C3) 



These expressions, coupled with (B7) and the fact that (f>o = (/)obs, fully define the emission coordi- 
nates. 
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For emission from within the cone of open field lines, several of these angles are small. If the 
observer is to see this cone, the difference between a and i must be less than (or on the order of) 
the colatitude of the last closed field line, ^c(^) = \/f/RL- Similarly, the magnetic colatitude of 
the emission, and the rotational phase ^ = Vlt must be of this same order. 

The expression for tan^o (B7) is perhaps the most difficult formula. It has two roots, and the 
continuous physical solution switches between them when ^0^5 = 7r/2. Since the maximum of ^of,^ 
is a + z, this switch will occur somewhere in the pulse phase for all pulsars with a + i > 7r/2. If 
we are only interested in the region close to the pulse center, then this is not a problem, but for 
millisecond pulsars and other pulsars whose emission covers a substantial fraction of the period, 
this switch-over has to be kept in mind. 

Here, however, for simplicity we assume that Oohs is always less than it/2. In that case, an 
excellent approximation to (B7) is 

tan 00 ^\ — ( 1 - cos dobs - 7^ sin^ Oobs ) (C4) 



while for small Oohs^ 



2 

tan ^0 ~ sin ^0 ~ o ^obs (C5) 



Expressing sin Oghs in terms of 

cos 6 obs gives us, for small Oobsi 

sin^ob, = [2(l-cos^o6,)]^/2 ((^g) 

which is quite useful, since we frequently have terms in even powers of sin^Oj and we know cos^obs 
from Eq. (C2). 

The observer and emission azimuth are the same, ^0 = ^obs- Using the transformed y- 
coordinate, ub = rsinO coscj), we find 

cos a sin z cos — sin a cos z 

cos<^o = C7 

sm Oobs 

. , -sinisinOi 

sm(^o = —a (C8) 

sm Oobs 

By using the relations between sin^o and sin^o^sj these can be cast into various useful forms. 
For example, in order to simplify expressions (2) and (6), we used these relations to write 



sin Oq'ih - sin O^bs ~ q { ^ ^ /? + - sin a sin i sin J (C9) 

and 

2 2/1 \ 

sin 00 cos ^0 ~ o ^i^ ^ohs cos ^ohs ~ q ( ^™ ''^ ~ o *^ ^^"^ ^ ) (CIO) 

where (3 = i — a and we have taken cosOt 1 — (1/2) sin^ As mentioned before, this limits 
these expressions to apply only to the region close to the pole where Ogbs is small. 
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D. Field Aligned Current 



We expect a current to flow along the open field lines with intensity approximately equal to 
the Goldreich- Julian charge density moving at the speed of light, 



J = -^C("-B)B 



(Dl) 



where ( = J/Jqj is scale factor of order unity, reflecting our ignorance of the actual current. 
Since we will be examining only first-order effects, the results for the perturbing fields and the 
polarization shifts are simply proportional to C,. In the interest of conciseness, we set C = 1 for the 
rest of this section and reinsert it into the final result. 

In the case of a small polar cap. the variation in • B across the cap is also small, so we take 
ri • B cos a. In addition, if B = Bq + Bi, with Bi <^ Bq and V x Bq = 0, as is true for any 

superposition of vacuum multipoles, then V x Bi = 2r2*Bo/c. 

Now assume Bq is a point dipole and introduce spherical polar coordinates with the z axis along 
the dipole axis and the origin of coordinates at the point dipole. Assume axisymmetry with respect 
to the magnetic axis. Since the assumed current is axisymmetric, this assumption is equivalent to 
assuming the polar flow tube has a circular cross section. 



The equations for the field components are 

1 d 



r sin 9 89 



{Bi^sm9) = - 



cos a cos 9, 



15 2n,fi . 

■-TrirBi^) = —cos asm 9, 

r or cr^ 



Id 1 dBir 

-^{rBie) — 

r or r o9 



0. 



Here Bor = {2/1/ r^) cos 9, Bqq = i/J'/r^) sin 9. The solution is 

Blr = BiQ = 0, 



B 



2/i / r 

— - cos a — 



sm( 



(D2) 
(D3) 
(D4) 

(D5) 
(D6) 



where Rl = c/fi*. One readily finds this Bi^p satisfies both of the non-trivial equations. 

This gives a perturbation to the tangent vector of 

r sin 6 

0, 0, -— cosa— — 
Kl iVi 

in (r, 9, (f)) polar components. 

Following the above method, and expanding in powers of sm9, the updated normal is 



(D7) 



ni 



r 1 
0, 0, -4— cos a(l - - sin^ 9) 
Rl 4 



(D8) 
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which results in a change in the binormal of 



bi = — — cos a 



-sine, 4(1 --s\n^e), 



(D9) 



The emission point changes by 



ri = -[Vto] ^ • ti = — cos a 



changing the observed binormal by 



5h = V6o • ri = cos a 



0, 0, -iVitan6' 



2 2 
gSin^, -,0 



giving a total perturbation of 



Ab = bi + (5b = -r— cos a 
Rl 



yielding a polarization angle shift of 



10 r J 



5 . . 10 .2 

— - smt^, — sm I 

3 ' 3 



7 



cos a 1 — — sin 



3 Rl Jgj V 40 
where we have re-inserted the dependence on the current strength. 



(DIO) 



(Dll) 



(D12) 



(D13) 



E. Return Current 

To model a return current, we assume that each field line has a fixed flux of current on it, 
proportional to the local GJ density. 



If the perturbing current density has the form 

j = _^ncosaB/(^) 
27r ■'^ r ^ 

for / an arbitrary function, then the perturbing field is 

2ni2 



B 



■ cos 



a^- / du' fiu') 



(El) 



(E2) 



where u = sin^ 6/r labels an individual field line. This is equivalent to using the magnetic flux as 
the integration variable. 



Dividing by the magnitude of the held gives 

Bi cos a r r 



\Bn 



Ni Rl sine 



r du'f{u') 
Jo 



(E3) 
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with Ni as defined in (B3). 



If f{u) = 1, we get the normal perturbing field. If we include a return current which is a 
multiple A of the background GJ current in a small ring past sin^ 6/r = ui, namely 



we get a perturbation of 



ti = — cos a < 



1, if «! > u > 0; 
-A, if (1 + 1/A)ui >u> ui; 
[0, if « > (1 + 1/A)ni; 



— sin 9, if ui > u > 

-^ii^^ + Asin^, if (1 + 1/A)ni > u > ui 



lo, 



if u> (1 + 1/A)ui 



(E4) 



(E5) 



This generates the normal polarization angle shift for u < ui, no shift for u > {1 + 1/X)ui, 
and a shift of 



10 r ^ fl + X uir 

Alp = — — Acosa — - 

6 Kl \ sm^ 6 

for (1 + l/A)'Ui > u> ui. If A S> 1, then this reduces to 



1 



(E6) 



10 r / 7 \ 

= "T ~ " 40 ''''' V ^^^^ 

where = X{u — ui) / ui ranges from to 1 through the return current layer. 

This simply eliminates the perturbation over the range from (1 + 1/A)tti > u > ui, which is 
what we expected. If the return current layer is thin, A » 1, this transition is abrupt. This would 
appear in the data as a constant negative shift outside of the region of current flow, as there would 
no longer be any current to offset the aberration, as shown in Figure 5. 



F. Aberration 

The aberrational perturbation arises because the velocity of particles streaming along field lines 
has components both along the rotation direction and along the field line. Following Blaskiewicz, 
Cordes, & Wasserman (1991), the velocity of the particles becomes 

v = c^l|B + J2xr, (Fl) 

where is a free parameter, reflecting the fraction of particle velocity along the field. Assuming 
the particles are moving at c and normalizing gives a change in the tangent vector of 

ft X r /Oxr 

to to (F2) 



-20- 



This change in the tangent vector is of the same form as the current aberration treated pre- 
viously, i.e. first-order in r/Ri and perpendicular to to- As such we can apply the same general 
formalism to calculate the expected polarization angle shift. Additional first-order changes to 
the magnetic field generate second-order changes in ti, so we may safely treat this perturbation 
independently. 



In polar coordinates, 

Cl X r = [0, sin a sin cf), sin a cos 9 cos cf) + cos a sin 6] 



(F3) 



Subtracting out the component along to gives 



r 



— ^N^'^ sin a sin cos sin (/) 

N^'^ sin a cos^ 9 sin <p 
sin a cos 9 cos d> + cos a sin 9 



(F4) 



Aberrational effects change the curvature of particle orbits in two ways. Not only does the 
instantaneous velocity of the particles change by f2 x r, but the underlying magnetic dipole is 
rotating as well. The adds a dto / dt term to the curvature vector perturbation k\ from ( A4) , giving 



dto 
dt 



+ (to • V)ti + (ti • V)to 



(F5) 



The normal vector is then found from (A5) to be 



ni 



r 



8 sin 
3 sin 



— sin a sin (l — | sin^ 9) 
— \ sin a sin 9 sin 
^(l-|sin2^)+4cosa(l-|sin2^) 



+ 0(sin3 9) 



producing a binormal change of 
bi 



r 



I sin a cos </> (l — ^ sin^ 0) -|- cos a sin 9 
l^^(l-i-n2^)-4cosa(l-|sin^^) 




+ 0(sin^ 



(F6) 



(F7) 



Solving for the emission point gives 



r 



— 3 Sin a sm ( 



■ sin2 9) 



— 3 sm a cos 



(l — I sin^ ^) — i cos a sin 9 



and a change in binormal of 



5h = 



Rl 



sm a cos ( 

2 sin a cos c, 

3 sin 5 



(1 



■ sm 



-9)-\ 

2 /i\ I 2 



+ 0(sin3 9) 



(1 



sin' 



cos a sm t 



cos a 



+ 0(sin=^ 9) 



(F8) 



(F9) 
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giving a total change in the binormal of 



35 
24 



sin^ ^) + t COS a sin ( 



sinacos(?!) (l — 

o sing cos (/> /-■ 47 „• 2 a\ 10 „^ ^ 33 2 a\ 

sing 11-24 Sl"^ - T COS « (l - 40 S™ 





+ 0(sin=^ e) 



and a shift in the polarization angle of 
r sin a cos (j!) /. 11 



Rj 



sin0 



sin^ I 



(FIO) 



10 r / 7 \ 

— — cosa 1-— sin^i9 +0(sin^6'). (Fll) 
6 Hl V 10 J 



G. Aberrational Phase Shift 

Qualitatively, the changes to the polarization angle caused by aberration appear to simply 
shift the phase of the sweep. BCW directly converted their calculated perturbation into a phase 
shift, but here we want to be slightly more precise and preserve the higher-order terms (in sin0). 
A portion of the aberrational change in the binormal Ab (FIO) corresponds to a phase shift, but 
part of it does not. In order to separate these two parts, we need to calculate how the binormal at 
the emission point changes with time. 

This will let us decompose the perturbation into Ab = A$(5b/5$) + Ab', which corresponds 
to a perturbed polarization sweep of -0 = V-'o(^ + ^^) + ^V"') where "00 (^) is the standard sweep, 
A-i/i' is calculated from Ab' as before, and $ = fii. 

As the star rotates, the zeroth-order binormal (as evaluated in magnetic coordinates) changes 
both because of the apparent motion of the observer across the sky and because of the change in 
the rotating coordinates themselves. The change in the binormal is 

dh dU ,idei db^^^ dlP'^ 

-r = -r^i + b -r^ = —r \ 1— (Gl) 

dt dt dt dt dt ^ ^ 

The first term depends on the motion of the emission point in magnetic coordinates. In 
these coordinates, the zeroth-order binormal only depends on the azimuth of the emission point, 
00 = (t'obs- From (C3), this changes in time as 

= — f2(cosa; + sinacos(?!)ocot^o6s). (G2) 

Solving for 9obs in terms of gives 

2-tan2 6'o 2 1 - sin^ 6*0 . 

cot Oobs = n ~ o + O sin3 9) G3 

Stanfc'o 3 sin^o 

The change in binormal as seen in the rotating frame is then 
"(1) 

db db d<po #0 

-dT = 5^^ = ^-""^°'-^°^^'^'']^ 
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= n 



2 sin a cos ( 



■ sin a COS (/)o (l — sin^ ^o) + cos a sin 9i 







3 singo (1 - i ^i^' ^o) + COS a (1 - I sin2 ^o) 




+ O{sm^eo) (G5) 



where we have expanded up to third order in sin^o- 



The second term depends on the rotation of the coordinate system. Defining Ry{a),Rz{b) to 
be respectively rotations around the y, z axes by angles a, b, this becomes 



dt 



(G6) 



or 



db 



(2) 



dt 



O, [— sin a cos 9 cos cp — cos a sin 9, — cos a cos 9 + sin a sin 6 cos ^, 0] 



(G7) 



Adding these two and expanding to third order in sin 9 gives 



*=0 

dt 



- - sin a cos 6 I 1 — - sin 



2 sin a cos ( 
' 3 sin^ 



-, 



+ 0{sm^9) 



(08) 



The shift in b due to aberration, (FIO), then corresponds to a phase shift of 

r 

A$ = -3 — 
Rl 

Once this phase shift has been subtracted out, we are left with a remainder 



Ab' = — 

Rl 



t cos a sin — || sin a cos cj) sin^ 9 



10 



cos 



a(l 



33 
40 



sm' 



''9) + 



sm a cos © sm ( 



+ 0(sin3 9) 



which corresponds to a polarization shift of 



(09) 



(OlO) 



AV'' = - 



10 r 
'Jr'l 



cos a [ 1 sin 

10 



47 r 

+ — — sinQ;cos0sin6' + O(sin^6') (Gil) 
12 Rl 
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